Ensemble of ultra-high intensity attosecond pulses from laser-plasma interaction. 
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The efficient generation of intense X-rays and 7-radiation is studied. The scheme is based on 
the relativistic mirror concept, i.e., a flying thin plasma slab interacts with a counterpropagating 
laser pulse, reflecting part of it in the form of an intense ultra-short electromagnetic pulse having 
an up-shifted frequency. In the proposed scheme a series of relativistic mirrors is generated in the 
interaction of the intense laser with a thin foil target as the pulse tears off and accelerates thin 
electron layers. A counterpropagating pulse is reflected by these flying layers in the form of an 
ensemble of ultra-short pulses resulting in a significant energy gain of the reflected radiation due to 
the momentum transfer from flying layers. 
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I. INTRODUCTION. 

The development of sources of intense ultra-short electromagnetic (EM) pulses, X-rays, and even 7-rays is an im- 
portant potential application of intense laser matter interactions Such applications vary from single molecule 
imaging to radiography of dense targets and medicine diagnostics, and provide the opportunity to study the funda- 
mental effects not available before, such as the effects of nonlinear Quantum Electrodynamics (QED). One of the most 
promising ways to generate such sources is the use of a reflection of EM radiation from a flying relativistic mirror. 
This was first studied by Einstein in [2] as an example of Lorentz transformations. The radiation frequency up-shift 
is proportional to the square of the mirror Lorentz factor, making the scheme very attractive for the generation of 
high frequency pulses. 

The principal idea of the relativistic plasma mirror has existed for a long time [3|]. Recently several ways to create 
such mirrors have been proposed. One way is to use the plasma waves in the wakefield of a high intensity pulse as 
it travels through low density plasma in the wave breaking regime [4]. The incoming light is reflected by these waves 
in the form of a compressed pulse with an up-shifted frequency (demonstrated in the experiment reported in Ref. 
(H|). Moreover additional intensification comes from the parabolic form of the wake wave [H, 01, which focuses the 
reflected radiation. The role of the counterpropagating electromagnetic pulse can be provided by different nonlinear 
structures in plasma, left in the wake of another pulse, such as solitons, electron vortices, or a wakefield. In this case 
the breaking plasma wave will reflect a part of their electromagnetic energy in the form of a single-cycle ultra-short 
high frequency pulse [8|]. Another potential method is the interaction of intense linearly polarized electromagnetic 
pulses with solid density plasma, where either sliding [§] or oscillating mirrors 0, Ell EH E2] can be formed. The 
part of the incoming pulse that is reflected carries high order harmonics due to the oscillation of the reflecting surface 
either in the transverse or longitudinal direction [H, Ell EH, E2| • 

It was proposed recently that relativistic mirrors can be formed in the regimes of laser-thin foil interaction previously 
considered in regard with the ion acceleration: in the first case the laser pulse is intense enough to separate the electrons 
from the ion core, so that the electrons will move in front of the pulse, forming a dense relativistic electron layer. The 
counterpropagating laser pulse upon reflection from such a mirror will undergo compression and frequency up-shift 
(l3l . EH- In the second case a double-sided relativistic mirror is formed. One side is used for the energy transfer 
from the laser pulse that accelerates the foil. While the other side can be utilized to reflect the incoming radiation to 
produce ultra-short high frequency pulses [151 ]. 

In this letter we propose a more realistic mechanism of generation of ultra-short EM pulses in laser-solid density 
target interaction, then previously discussed in the literature fl3l EH- Sice it represents a more general scenario 
of ultra-short pulse generation in laser-thin foil interaction. In the proposed scheme the ultra-short EM pulses are 
created in the course of a counterpropagating laser pulse interaction with a series of flying electron layers, i.e. the 
Relativistic Multilayer Reflection (RMR) mechanism. Such layers are produced when an intense laser pulse interacts 
with a thin solid density target and extracts and accelerates thin electron layers [H, E3, El, EH • The high density and 
relativistic velocity of these electron layers make it possible that such structures will reflect the incoming radiation, 
acting as flying mirrors. Since the incident pulse experiences the Relativistic Multilayer Reflection, such interaction 
will result in the generation of an ensemble of ultra-short pulses. Moreover this scheme can lead to a significant energy 
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gain by reflected radiation due to the momentum transfer from flying layers. In what follows we first present the 
results of 2D PIC simulation of the proposed mechanism and then a ID analytical model is employed to analyze the 
results. 



II. THE RESULTS OF 2D PIC SIMULATIONS. 



In our numerical model with the 2D PIC code REMP - relativistic electromagnetic particle - mesh code based on 
the particle-in-cell method [20] the generation of ultra short pulses by reflection from the series of laser accelerated 
dense electron layers is studied in high-intensity laser interaction with ultra-thin targets. The targets are composed 
of fully ionized carbon C +6 with an electron density of 400n cr . The grid mesh size is A/200, space and time scales are 
given in units of A and 27r/cc;, respectively, the simulation box size is 15A x 12. 5A, where A and uo are high-intensity 
laser wavelength and frequency respectively. The number of particles per cell is 225. The 1.6 PW laser pulse which 
generates flying electron layers is introduced at the left boundary and propagating along x axis from left to right. The 
pulse is linearly polarized along the y axis (P-polarization), tightly focused (f/D = 1), and has Gaussian transverse 
and longitudinal profiles. The duration of the pulse is 30 fs. The counterpropagating laser, which is reflected from 
the relativistic mirrors, is introduced at the right boundary and propagates from right to left along the x axis. It is 
polarized along the z axis (S-polarization) in order to distinguish between the radiation generated by the accelerating 
pulse and the reflected one. It has ao = 1, Ao = 4A, and duration of 15 fs. 

Below we present the results of 2D PIC simulations for the cases of a mass limited target and a thin foil (see Fig. 
1). In the first case the counterpropagating laser pulse is reflected from the flying electron layers accelerated from a 
disk with diameter 1A and thickness 0.1 A placed at x = 6. OA, i.e. before the focus of the accelerating pulse, which is 
focused at x = 8. OA (Figs. la-c). Such configuration is chosen in order to show the formation of flying electron layers 
without a transverse flow of electrons towards the irradiated spot and obtain a more regular group of flying electron 
layers. In the second case the pulse interacts with a 0.1 A thick foil placed at the focus (x = 6 A) of the accelerating 
pulse. The electron density distributions at t = 20 are shown for both cases in Figs, la and Id. The thin electron 
layers that act as flying relativistic mirrors can clearly be seen. The density of these layers vary from 5 to 15n cr for 
the mass limited targets and from 10 to 20n cr for the thin foil. The duration of these bunches is about 70 as. 
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FIG. 1: (color on-line) The reflection of laser pulse by accelerated electron layers in the case of a mass limited target (upper 
row: a, b, and c) and a thin foil (lower row: d, e, and f). The electron density distribution a) and d) at t = 18; the distribution 
of counterpropagating pulse electric field after reflection, b) and e) at t = 18, the field is measured in units of m e cco/e; the 
distribution of electrons in (p x ,x) phase plane at t = 12, c) and f). 
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III. ID MODEL OF RELATIVISTIC MIRROR GENERATION AND RADIATION REFLECTION BY IT. 



Let us estimate the properties of the flying electron layer using a simple ID model. As we have seen from the results 
of 2D PIC simulations the electron layers are generated by the consecutive maxima and minima of the laser electric 
field. The number of electrons per flying layer can be obtained from the requirement that the Coulomb attraction 
force should be greater than the Lorentz force exerted by the EM field on these electrons. Since all the extracted 
electrons almost instantly become relativistic, v e = c, the Lorentz force is proportional to laser electric field. Or in 
terms of fields, the laser field should be greater than the charge separation field, E cs = 7re5n e l, which is generated after 
the extraction of electrons (I is the foil thickness, Sn e is the density of evacuated electrons). This gives the condition 
a > 7r5n e l/n cr \, first introduced in Ref. [11] as a condition for thin foil transparency. Here a is the amplitude of a 
vector potential in the corresponding maxima or minima of the field, n cr is the critical plasma density, and A in the 
laser wavelength. Here we should note that the density of evacuated electrons is Sn e if Sn e < n e , otherwise all the 
electrons are evacuated and Sn e = n e . 

However, as it was pointed out in [19], this condition does not take into account the fact that each electron layer, 
escaping the attraction of the ion core, increases the charge separation field that should be compensated for by the 
laser pulse field. Let us approximate the field of the laser as ao exp[— t 2 /r 2 ] cos[27r£/T], where r is the half of the 
duration of the pulse and T is the period of the EM wave. The maxima and minima of such a wave are at tj = Tj/2, 
j = 0,±1,±2, .... Then for some field maximum ctj the charge separation field that already exists is determined by 
aj+i. Then the number of electrons evacuated is A7VJ = \R 2 n cr (a 3 ; — a^+i), here R is the radius of the irradiated area. 
The total number of evacuated electrons will be A7V e = J2j ~ Ai? 2 n cr ao, i.e. is determined by the maximum of 
the vector potential only. Here we also assume that A7V e < 7V e , otherwise A7V e = N e . 

Let us estimate the duration of an electron bunch. The extraction of electrons begins only when a > a J+ i and stops 
at the top of the current laser cycle a = dj. Then the time interval £ that determines the duration of the bunch can 
be found by solving the equation a(tj + £j) = cij+i* 
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The duration is minimal at the maximum of the pulse (j = 0) £,o/T = 0.11(T/r) and increases with the increase of j. 
For a laser pulse with r = 5T and T = 3 fs the electron bunch has an attosecond duration: £o = 60 as, which is in an 
agreement with the results of 2D PIC simulations as well as the number of electrons per bunch. 

Such thin flying electron layers can reflect the conterpropagating radiation in the form of short pulses with up- 
shifted frequency. The generation of ensembles of such EM pulses through the RMR mechanism is demonstrated in 
Figs lb and le. In order to determine the properties of flying electron layers and determine Lorentz factors of mirrors 
the distributions of electrons in (p x ,x) phase plane are shown in Figs, lc and If for t = 12. The formation of flying 
relativistic mirrors with 7^5 can clearly be seen. 

In order to characterize the generation of ultra-short EM pulses we show the results of spatial Fourier analysis of 
reflected (curve 1) and incident (curve 2) radiation in Fig. 2 for the case of the mass limited target. The calculations 
are preformed for the field along the line ?/ = 0.6Aatt = 15, when the maximum frequency radiation can be resolved 
on our grid. As the interaction evolves we expect further frequency up-shifting based on our theoretical model. It can 
be clearly seen that high-frequency radiation is generated as a result of counterpropagating pulse reflection from flying 
electron layers. The spiky structure of the spectrum can be a consequence of the regular distribution of reflecting 
electron layers in space (see Fig. la). 

Let us estimate the reflection coefficient of the flying electron layer to determine the efficiency of short pulse 
generation mechanism. In order to do so we perform Lorentz transformation to the reference frame co-moving 
with the electron layer and use the results of Ref. [11], where the interaction of an EM wave with a thin foil was 
considered and the reflection (p = €j/(i + e^)) and transmission (r = + e^)) coefficients were obtained. Here 
ej = 6q/(1 + /?)7 = 6q/27 (since /3 ~ 1) is the parameter governing the transparency of the foil in the co-moving with 
the foil frame and e J = nAn^j /n cr T is the transparency parameter of the electron layer in the laboratory frame 
Since it is essential for our mechanism that the flying electron layers remain unperturbed by counterpropagating pulse 
we take the intensity of such pulse to be relatively small. 

Let us estimate the intensification of the laser pulse reflected by such a mirror. For the incident laser intensity, Iq, 
the reflected pulse in the laboratory frame will have the up-shifted frequency by a factor of 4j 2 and the increased 
intensity 
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FIG. 2: (color on-line) The results of spatial Fourier analysis of the reflected (l,red) and incident (2, black) pulses in the case 
of a mass limited target. 
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For 7 <C 6q the reflection coefficient tends to unity and reflected intensity end energy are determined by the Lorentz 
factor alone. In the second case, 7 ^> ej, the foil moves so fast that it becomes increasingly transparent for incoming 
radiation and the energy of the reflected pulse is limited by the transparency parameter. Thus the efficiency of the 



light intensification is determined by 7 2 x min{7 2 , e 3 } for intensity and min{7 2 , e J } for energy. This means that the 
maximum intensity increase is determined by the properties of the accelerated foil. In order to illustrate the results 
of the theoretical model we chose e J ~ 10 and 7 ~ 10 as parameters to estimate the reflected pulse intensity. These 
numbers come from the results of 2D PIC simulations. The frequency up-shift in this case is 4 x 10 2 and for T = 3 fs 
the period of the reflected radiation will be T r = 7.5 as. Then the intensity is I ~ 3 x 10 22 W/cm 2 . 

As we have shown the energy of the reflected from the flying electron layer EM pulse is limited by e J 2 , if we 
consider the dependence on 7. However there is a way to increase the energy of reflected radiation by utilizing the 
fact that the accelerating laser pulse extracts multiple electron layers from the target. Moreover this type of reflection 
is realized in the results of 2D PIC simulations presented above (see Figs la and Id). The counterpropagating pulse 
experience relativistic multilayer reflection giving rise to an ensemble of ultra-short pulses. While the number of 
photons reflected at each layer, AN^ = |p| 2 A^ (N® is the number of photons in the counterpropagating pulse) can 
be small, the multiple layers will reflect almost all incoming photons, — > N®. If the number of layers is /, then 
the total number of reflected photons [21] is 



n: 



=^i P i 2 x;(i-H 2 ) i =^[i-(i-ipiv]. 



Here we assumed for simplicity that all electron layers are the same. Then in terms of energy 

Sr = 4 7 2 [1 - (1 - \p\ 2 y] So- 



(4) 



(5) 



If \p\ 2 f <C 1 then S r = 47 2 |p| 2 /f . In the limit / — > 00 the reflected back energy equals to S r = 47 2 £ and = N®, 
i.e. all the counterparpagating radiation is reflected. The energy gain is due to the momentum transfer from the 
flying electron layers to the reflected radiation. For the parameters of the electron layers obtained in simulations (the 
average density and duration of about 10n cr and 60 as, the number of layers is equal to 4) the reflection coefficient 
is 0.05 and the energy gain is 1.5, according to ([5]). We can estimate the actual energy gain in the results of 2D PIC 
simulations from the spectra of incident and reflected radiation by integrating \E z {k x )\ 2 over dk x . This gives a 1.07 
energy gain. This clearly shows that the proposed mechanism of RMR generation of ultra-short pulses leads to an 
energy gain of reflected radiation due to the momentum transfer from flying mirrors. 
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IV. TOWARDS THE SCHWINGER FIELD. 



Further intensification can possibly be achieved by focusing by some external means the reflected pulse into a 
diffraction limited spot. For a single pulse it will lead to 



where D is the reflected pulse width before focusing and A r = Xo/4j 2 . Then for the parameters of the flying layer 
used above, D = 3Ao and io ~ 10 18 W/cm 2 the resulting intensity will be of the order of the intensity characteristic 
for the effects of nonlinear QED, i.e. Schwinger intensity, Is ~ 10 29 W/cm 2 , [22j]. At this intensity the probability 
of one of the most profound processes of nonlinear QED, the e + e~ pair production in vacuum by strong EM field, 
becomes optimal. We should note here that the plane EM wave does not produce pairs in vacuum [22|, because in 
this case both field invariants, T = E 2 — H 2 , Q — EH, are equal to zero, which is not the case for the focused pulse 
(23) . Thus the focusing of high frequency pulses to ultra-high intensity will provide a unique opportunity to study 
both the e + e~ pair production by a focused EM pulse and its dependence on frequency of the pulse [24|]. 



In this letter we considered a new way to generate ultra bright high intensity X-rays and gamma-rays by reflecting 
EM pulse from the relativistic mirror. In the proposed scheme the role of the flying mirror is taken by laser accelerated 
electron layers, which are formed in the process of the intense laser pulse interaction with thin solid density targets. The 
reflected pulses have an up-shifted frequency and increased intensity. The reflected pulse intensification is determined 
by a combination of two parameters: the Lorentz factor of the flying electron layer, 7, and its transparency parameter, 
[HI- It is proportional to 7 2 x min{7 2 , (e J ) 2 } and for given density of the plasma slab the energy of the reflected 
pulse multiplication factor can not exceed (e J ) 2 . It is due to the fact that for 7 ^> e J the fast moving plasma slab 
becomes increasingly transparent for the incoming radiation and the amplitude of the reflected pulse drops, limiting 
the energy gain. Further intensification of the reflected light can possibly be achieved by its focusing into a diffraction 
limited spot that will bring the resulting peak intensity well into the domain of nonlinear QED with laser systems, 
which are presently available. 

We showed that there is a way to increase the energy transfer from the accelerated layers to the reflected radiation 
by utilizing the relativistic multilayer reflection. The counterpropagating EM pulse interacts with multiple flying 
electron layers, producing an ensemble of ultra-short pulses with an energy scaling of 47 2 . This fact leads to the 
conclusion that this mechanism of ultra-short pulse generation provides a highly efficient way of transforming the 
incoming pulse into high frequency radiation through the relativistic multilayer reflection. 

This work was supported by the National Science Foundation through the Frontiers in Optical and Coherent 
Ultrafast Science Center at the University of Michigan and Russian Foundation for Basic Research. 
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